Self-similar solutions with fat tails for a coagulation 
equation with diagonal kernel 
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Abstract 

We consider self-similar solutions of Smoluchowski's coagulation equation 
with a diagonal kernel of homogeneity 7 < 1. We show that there exists a 
family of second-kind self-similar solutions with power-law behavior x-^^+p'^ as 
x — >■ 00 with p S (7, 1). To our knowledge this is the first example of a non- 
solvable kernel for which the existence of such a family has been established. 

1 Introduction 

Smoluchowski's coagulation equation provides a mean-field description of binary 
coalescence of clusters. If denotes the size of a cluster and fi(,,t) the corresponding 
number density at time t then the equation is 



f{i,t) = \l drjK{^-r],r])f{r,,t)f{^-r],t)-f{^,t) I drj K{^,r])f{r,,t) , (1) 



dt - JO JO 

where K{(^, rj) is a kernel that describes the rate of the coalescence process. 

Here we consider a specific diagonal kernel of homogeneity 7 < 1, given by -fC(^, r/) = 

5(^^_fj^^^~^'^ , that reduces ([I]) to 



dt' ' 4 V2' / ' V2' 
In the following we study self-similar solutions of ([2]). Such solutions are of the form 

^(l+(l+7)/3) J ^ 



fi^,t) = t'''^''^'^'^g[j^) (3) 



for some positive /3, where g satisfies, with x = S,/t^, that 

1 /x\ 1+7 /x^ 



;i + (l + 7)/3)<7-/3x5'(:c) = i(|) V(f)-xi+V(x). (4) 
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If one looks for solutions with conserved mass, then (3 is uniquely determined by 
(3 = /3^, := 1/(1—7). For further reference we also note that we can integrate the 
equation in ([ID to obtain 

^x^gix) = r s^+^g\s) ds + (l-7)(/3 - f sg{s) ds . (5) 

Jx/2 Jo 

Here we assumed implicitly that xg{x) and x'^'^'^ g^ {x) are integrable at zero and that 
\mix^Qx'^g{x) = 0. As we will see below (cf. @), these properties will be satisfied 
by the solutions we are going to consider. Notice also that we have the well-known 
power-law solution 

^^-(1+7)^ with ^ := 2^-1 < 1 . (6) 

In [1] a mass-conserving solution of ([5]), that is a solution for /3 = is constructed 
that is decaying exponentially fast and satisfies 

^(x) =x~(i+^)(^ -cx'^/^i'^) +o(x'^/(i~^))) asx^O, (7) 
V 1 — 6 J 

where /i > satisfies a certain transcendental equation. The constant c > is not 
determined due to an invariance of dH) under the rescaling g{x) 1— )• a^^^g(ax) for 
any a > 0. In the case of mass-preserving solutions the constant can be fixed by 
normalizing the mass of the solution. As is pointed out in [1], the solution is unique 
in the class of functions satisfying ([7]), but uniqueness in general is not known. 
In [1] the question is raised whether solutions with algebraic decay, others from 
the one in ([6]), exist in analogy to the ones that have been found in p] for the 
constant and additive kernel. More precisely, for example for the constant kernel, it 
is established in [2] that there exists a family of self-similar solutions with infinite 
mass and the decay behavior x~^^^^^ for all p G (0,1). Furthermore, it is shown 
that a solution of the coagulation equation converges to the self-similar solution 
with decay behavior x~^^^^^ if and only if the mass-distribution of the initial data is 
regularly varying with exponent p. In this note we prove for the diagonal kernel the 
existence of a corresponding family of self-similar solutions with infinite mass and 
asymptotic behavior x~^^~^p^ as x ^ 00 with p G (7,1)- Notice, that this includes 
solutions that are increasing asx— >ooif7<— 1. Our proof is simple and exploits 
strong monotonicity properties of a suitably rescaled version of the equation for the 
self-similar solution. We presently do not know, however, how to characterize the 
domains of attraction of these self-similar solutions. The analysis in [2] relies on the 
fact that the Laplace transform of the equation satisfies a simple ODE, a method 
that is not applicable in the present situation. 
Our main result is the following. 
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Theorem 1.1. Let 7 < 1 and fx be the unique positive solution of 

1 + /3;u _ 1 - 2^-1-^ 
2 ~ 1 - 2T-1 ■ 

Then there exists for any (3 > a solution g of ^ such that 



(8) 



g{x) = - CX^/(1~T) + o(x^/(l-^))) (9) 

as X ^ with a positive constant c Furthermore, x^^^^^"^ g{x) is monotonically 
decreasing and satisfies 

~ asx^oo (10) 

for some positive constant d. 

As explained above, the constants c and d in Theorem 11.11 are not determined due 
to the invariance of the equation under appropriate rescahng. 



2 Proof 

Our proof proceeds similarly to the one in [1] for the mass-conserving solutions. 
First, to scale out the singular behavior as x — t- 0, we introduce h{x) = g{x)x^'^'^ 
such that h solves 

or, due to ([5]), 

I3x^-^h{x) 

Notice, that the power-law solution ^ corresponds to the constant solution h = 
1/(1—^). It is also clear that any solution of (jlip for which \\mx^Qh{x) exists, that 
this limit must equal 1/(1—0). We are now looking for solutions that bifurcate from 
this constant at x — )• 0. 

In order to identify the next order behavior, we make the ansatz h{x) = 1/(1—0) + 
+ o(x'^) as X ^ 0. Plugging this into p^ . recalling that /3* = 1/(1—7) S'^d 
rearranging we find that fi must indeed satisfy ([8]). If we denote by F{fi) = (1 — 
2'^~^~^^) / il—O) we see that -F(O) = 1 > 1/2. On the other hand, F is increasing and 
lim^_j.oo Fil^) = 1/(1~^)- Hence, there must be a unique positive solution of ([8]). 
Next, we introduce the function j{x) via 



I3xh'{x) - h{x) = 6h^ (I) - h^{x) 



(11) 



/ s-"'h\s)ds + {l--/){/3 - /3,) / s-"'h{s)ds. (12) 

Jx/2 Jo 



= YZ^ + ^^(-c + j(2;)) , 



(13) 
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where c G M is a constant. Using equations ([8|) and ()12p we obtain that j satisfies 
1 



+ (l-7)(/3-/3,) / s-^+'^jis)ds)=:T[j]. 



In order to prove that a local solution of (I14p exists, we can proceed analogously to 

[1]. We only indicate the main steps here. 

We define for some e G {0, fi) and z > the space 

Q(z):={/eC[0,z];/(0)=0; := sup x~'\fix)\ < <x,} . 

It is clear that the operator T maps Ce{z) into itself. Next, we are going to show that 
T maps a ball in Ce{z) of a sufficiently small radius R into itself if z is sufficiently 
small. This follows from 

||T[j]|| < -mi— i + lljll — (1 _ 27-1-/^-2^) 

I ^2 , (l-7)(/3-/3.: 



1-7+2^ l-7+/i+e 
that implies 

Now we know by the definition of fi that 2F{fi + e) < 1 + f3{fi + e) and hence 
^ -(2F(^ + e) + (l-7)/3 - 1) < ^_(;,+e+l-7) = 1 . 



/3(l-7+/i+e) ' '"^ 1-7+^+e 

Thus, there exists a constant k = K(e) < 1 such that if ||j|| < R we find < 
kR + Cz'^{R'^ + !)• For sufficiently small z and an appropriately small R the right 
hand side is bounded by R. Similarly one can show that T is a contraction, we omit 
the details here. Hence, a local solution to ([HD exists, and thus also to p!T|) . Next, 
we choose c > 0, and claim that h is decreasing in a neighbourhood of zero. To see 
this, notice that it follows from ()14p that j'{x) exists for x > and that we have the 
estimate < (7-1^^^ + Cx^~^ for x G (0,z). This in turn implies that 

h'{x) = -c + + xV(^) < x''-^{- c^l + C^\j{x)\ + Cx>') . 



If z is sufficiently small, we find that h'{x) < for x G (0, z). We are going to show 
that as long as h exists and is positive this property is conserved. Indeed, assume 
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that there exists xq > such that /i'(xo) = 0. Then (jlip and the fact that h is 
decreasing for x < xq imply that 



= h{xof - h{xo) - ^^'(y)' < {l-0)h\xo) - h{xo) = h{xo){{l-e)h{xo) - l) . 

As long as h is positive, the right hand side is strictly negative, since h{xo) < 1/(1—6') 
and we obtain the desired contradiction. Moreover, equation (fT2|) implies for /? > /3* 
that h is positive whenever it exists. Hence, using standard results on ordinary 
differential equations, we obtain global existence of a solution h to (jlip which is 
strictly decreasing. Since equation (fTT]l has the only stationary points 1/(1—0) and 
0, it also follows that h{x) ^ as x — ?> oo. 

It remains to show that h{x) ~ dx^^^^ as x ^ oo from which (llOp follows. First, due 
to the invariance of equation (jlip under the transformation x — >■ ax for a > 0, we can 
assume without loss of generality that h{l) = 1/2. Since h satisfies f3xh'{x) + h{x) < 
h?'{x) we have by simple comparison that 

h{x) < ^-TTT for X > 1 . (15) 

We now introduce p{x) = x^^^h{x) that solves 

/3p'{x) = x-(i+i/'')(/(x) - 022/V(|)) . (16) 



The estimate (llSp in particular implies that p{x) < 1 for all x > 1 and thus (jl6p 
implies that P\p'{x)\ < 2x~(^"'"^/'') for all x > 2. Hence \p{x) — p{xo)\ < 2xq^^^ for 
any xq > 2 which implies that limx^ooP{x) exists. In order to complete the proof of 
Theorem 11.11 it remains to establish that this limit is strictly positive. To this end 
we note that (fT2|) implies 



X 



/3x^-^/i(x) > (l-7)(/3-/3,) / s-^'h{s)ds. (17) 







If we define ^ (x) := s '^h (s) ds then pT|) implies that /3x<I>' (x) 
— (1 — 7) (/3 — /3*) $ (x) > 0. Integrating this last inequality we obtain 



(l-7)(/?-/3.) ' ' 



<I> (x) ) > and thus 

$ (x) > $(1) = / s~'^/i(s) dx =: Co > 
^0 



_{i-7)Ca-A0 
X /3 



for all X > 1. Thus 



(l-7)(/3-/3*) 1 _1 

<I> (x) > CoX P = CqX ^x 
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for X > 1 and plugging this into (fT7|) we find h{x) > for all x > 1, that 

finishes the proof. 
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